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We determine the motions of the roots of the Bethe ansatz equation for the 
ground state in the XXZ spin chain under a varying twist angle. This is done by 
analytic as well as numerical study at a finite size system. In the attractive critical 
regime < A < 1, we reveal intriguing motions of strings due to the finite size 
corrections to the length of the strings: in the case of two-strings, the roots collide 
into the branch points perpendicularly to the imaginary axis, while in the case of 
three- strings, they fluctuate around the center of the string. These are successfully 
generalized to the case of n-string. These results are used to determine the final 
configuration of the momenta as well as that of the phase shift functions. We 
obtain these as well as the period and the Berry phase also in the regime A < — 1, 
establishing the continuity of the previous results at —1 < A < to this regime. 
We argue that the Berry phase can be used as a measure of the statistics of the 
quasiparticle ( or the bound state) involved in the process. 
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1. Introduction 



The one-dimensional spin one-half XXZ chain with the Hamiltonian 



H = 



1 N 



(i.i) 



provides a soluble many-body problem which has been intensively studied for more 
than two decades and is the subject of the present paper* Here a^ ,y,z are the Pauli 
matrices at site j and A denotes the parameter of anisotropy. Let a generic M-body 
state of the model be represented by 



N N M 

Xl = l Im=1 j=l 



1.2) 



where | f) is the state with all spins up and a x is the spin- lowering matrix at site 
x. The eigenfunction x( x ii ■ ■ ■ > x m) is of the form proposed by Bethe [1]: 



x(xi,...,x m \pi,---,Pm) 

X a ) 



IJ e(x b 

M>b>a>l 



El 





cxp 



M 



o=l 



^ X aPQa + ^ 0(pQb, PQa)t(Xb ~ X a ) 



M>b>a>l 



(1.3) 

where e(x) is the sign function, the summation is with respect to permutations Q 
of the momenta p and 6? is the two-particle phase shift. For two particles with mo- 
menta p and q, the explicit form of the phase shift is 9(p, q) = 2arctan 



Asin[(p-g)/2] 



cos[(p+g)/2]-Acos[(p- 9 )/2]_ 



The energy of this eigenstate is found to be E = — + 2 ^^(A — cospj). The 
second term is the sum of the single-particle energies Sj = 2(A — cospj). The 
operator M = — °f ) i s conserved and the case of half-filling corresponds 



* We find it impossible to exhaust all references on this subject. The partial list includes 
[1]-[12]. The developments in recent years include the relation to the study of finite size 
corrections as well as the approach based on the afhne U q (sl(2)), which are seen, for example, 
in [8], [9], [11], [12]. 
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to M = N/2. The complete reducibility of the problem to the two-body problem 
with the dynamical phase shift as well as the rigorous notion of the quasiparticle 
excitations is the hallmark of the integrable quantum systems whose mathemati- 
cal structure is captured by the Yang-Baxter relation [13] [14] and the attendant 
quantum inverse scattering formalism [15]. 

The boundary condition is a heart to this problem as this determines the 
spectra of the model and is referred to as the Bethe ansatz equation in this paper. 
Let us now impose, contrary to the standard periodic boundary condition, the 
twisted boundary condition 

X(xi, . . . , Xj + N, . . . , x M ) = e l ®x(xi, ...,Xj,..., x M ), (1.4) 

where $ is a twist angle.^ This leads to the system of transcendental equations for 
the set of momenta {pj} 

M 

N Pj + 0(Pj,Pk)=^Ij + ®, (1-5) 

fc=i(#) 

where j = 1, . . . , M. We denote by {Ij} the set of half-integers for even M or 
integers for odd M. 

In this paper, we study ground state properties of the system when the twist 
angle $ is varied. The change of an external parameter like $ in integrable systems 
will correspond to an adiabatic change in more general systems which do not possess 
an infinite number of conservation laws . The process we consider may be called 
adiabatic process in this sense. The major purpose of this paper is to examine the 
motions of the roots of eq. (1.5) under this process and to determine the momenta 
and the phase shift functions for the final configuration, given the ones for the 
initial configuration. 

f The twist angle plays an important role in projecting out some states to construct another 
model. See [10] [16] for example. 
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It is well-known that the system can be represented by the Jordan- Wigner 
fermions. The state with spin-up is regarded as an unoccupied state and the 
state with spin-down as an occupied state. The annihilation operator Cj and the 
creation operator cj are respectively defined by Cj = ^(<rj + ia y -) Y\k<j °t' Cj = 
\{<y X j — icTj) Ylk<j °t> which satisfy anticommutation relations {Cj, C\} = 8j^. The 
Hamiltonian (1.1) is rewritten in terms of lattice fermions as 

3 3 

The twisted boundary condition can be attributed to the twisting of the operators 
Cj = e i<s> C j+N , C] = e~^C] +N . The conserved operator M = J2f =1 CjC) 
measures the fermion number. The twisted boundary condition is related to the 
Aharonov-Bohm effect if we consider the M fermions carrying a charge (—q) on 
the ring threaded by the magnetic flux as each fermion goes around the ring, 
it picks up the Aharonov-Bohm phase $. 

The change of the twist angle $ is equivalent to varying the magnetic flux 
penetrating the ring which generates an electric field around the ring. The whole 
spectrum of Hamiltonian is periodic with respect to $ with period 2-7T. If we 
follow each individual energy level, this is not necessarily the case however. The 
period of the ground state, in particular, is a nontrivial quantity to be determined. 
In the case of no interaction A = 0, the period can be obtained by a classical 
consideration alone. All negative energy pseudoparticles are accelerated in the 
electric field E generated by the change of the flux A<3> through A$ = $(£/) — 
$(tp) = J* 1 dt <f dxqE . The change of the momentum is Apj = dtqE = ^ 
where N = § dx. Thus the period of the ground state is A$ = 2Nir. In this process 
all negative energy pseudoparticle move from the first Brillouin zone — n < p < n 
to the second Brillouin zone n < p < 3n. 

A completely different picture is known to be obtained as soon as the the 
interaction is turned on [17]. In [17], the regime — 1 < A < is considered. It is 
found that only one negative energy pseudoparticle moves from the first Brillouin 
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zone — 7r < p < n to the second Brillouin zone n < p < 3n. The period of the 
ground state is found to be A$ = An. During this process a level crossing occurs. 
The ground state becomes the excited state for the system with the original periodic 
boundary condition at $ = 2n and comes back to the original ground state at 
$ = An. The Berry phase attendant with this process has been considered in ref. 
[18] and is found to be n. 

In the regime < A < 1, more interesting events take place. It is pointed 
out in ref. [19] that, in the region cos(7r/n + 1) > A > cos(7r/n), a bound state 
of the n-string is formed during the process. Only these roots of the Bethe ansatz 
equation forming the n-string move from the first Brillouin zone —n < p < n to 
the second Brillouin zone n < p < 3n. The period of the ground state is again 
found to be A$ = An. 

The thrusts of the present paper consist of the following several points. In 
section 2, we consider the region < A < 1. In the segment cos(-7r/n + 1) > 
A > cos(7r/n) for n — 2, 3, • • •, we determine the motions of the roots of the Bethe 
ansatz eq. (1.5)in the rapidity plane*. This is done on the basis of our numerical 
and analytic study at a finite size system. Here we find intriguing motions of the 
n-string during the process. In order to state our results here more clearly, let us 
recall that the n-string is a group of rapidities which differ in their imaginary parts 
by 2i(n — ji) and center on the in line (ImA = n): 

= [A(»> + in] +i(n-n)(n + l- 2k) + 6™ 
k — 1, . . . , n 

where [\^ + in] is the center of the n-string and 5^ is a deviation due to the 
finiteness of the system from the string found in the infinite volume limit? This 
deviation is at most of order 1/N. We find, however, that it is this deviation which 

* In this paper, the words "root" and "rapidity" are used interchangeably, 
f In this paper, we deal with a natural formation of the n-string configuration. The Takahashi 
condition [6] will be automatically satisfied. 
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plays a vital role in determining the global structure of the motion of the roots 
under the process. To be more explicit, we find that, in the case of the two-string, 
the roots collide into the branch points of the momenta and those of the phase 
shift functions perpendicularly to the imaginary axis. In the case of the three- 
string, we find that the roots fluctuate around the center of the string due to the 
finite size correction 5^ mentioned above. The case of the n-string is found to 
a generalization of these: for n odd, the roots oscillate around the center of the 
string while, for n even, in addition to the fluctuation, the two roots in the middle 
collide into the branch points. This fluctuation of string state creates a winding 
around the cut of the momentum and the phase shift function and determines the 
final configurations of these. These are presented in section 2. 

In section 3, we consider the regime A < — 1 and determine the motion of the 
roots. From this we determine the period, the final configurations of the momenta 
and the phase shift functions and the Berry phase. These values are found to be 
the same as those derived in ref.[17], [18] in — 1 < A < 0. We, therefore, establish 
the continuity at A = — 1. 

In section 4, We compute the Berry phase in the regime < A < 1 and consider 
the physics associated with it. By piecing through the computation done in [18] and 
ours, we are led to propose that the Berry phase provides a measure of the statistics 
of the quasiparticle (or the bound state) involved in the process. We consider 
statistics of particles as the Berry phase arising through p-space monodromy. 

In Appendix A, we summarize our analytic study of the finite size system at 
M = 2, N = 4. In Appendix B and Appendix C, we present proofs of technical 
lemmas necessary in §2 — 1 and §4 respectively. 
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2. Motions of String Solutions under a Varying 
Twist Angle in the Regime < A < 1 

In the regime — 1 < A < 1, the momentum p, the phase shift 8 and the total 
energy E are parametrized by rapidities Xf 

p(X) = — im 

0(Ai,A 2 ) =i\n 

E — -2sm z /x 

2 C ° S ^ coshAj — cos/i ' 

where A = —cos /i. Let us imagine a closed circuit in the A plane consisting of the 
real axis and the in line connected at plus and minus infinity. The momentum p(X) 
increases monotonically as a function of A from — (n — p) at A = — oo to {n — p) 
at A = oo. ( See Figure 2.1). If p increases further, A moves backward along the 
in line from +oo + in. The momentum p(X) varies from {n — p) at A = oo + in to 
2n — (n — p) at A = — oo + in. When the rapidity goes around this closed circuit, 
the momentum p increases by 2n, since the path encircles a logarithmic branch cut 
as a function of A. We note that the single-particle energy 

— 2sin 2 /i 
3 coshAj — cos/i 

is negative for Im Xj = and positive for Im Xj = n. 

The behavior of the phase shift function 6(X) depends on the sign of A. ( 
See Figure 2.1.) When A < (0 < p < |), the phase shift 0(A) decreases from 
(n — 2p) at A = — oo to — (n — 2p) at A = oo. The phase shift further decreases 
from — {n — 2p) at A = oo + in to — 2n + {n — 2p) at A = — oo + in. As the 

8 



sinh^(A — ip) 
sinhi(A + ip) 



(2.1) 



sinh^(Ai — A2 — 2ip) 
sinh^(Ai — A2 + 2ip) 



(2.4) 



relative rapidity goes around the circuit, the phase shift decreases by 2n. When 
A > < // < 7r), the phase shift 0(A) increases from — (2/x — 7r) at A = — oo to 
(2/x — 7r) at A = oo. The phase shift further increases from (2/x — it) at A = oo + in 
to 27r — (2/x — 7r) at A = — oo + «7r. As the relative rapidity goes around this circuit, 
the phase shift increases by 2n. This difference plays a central role in determining 
the behavior of the rapidities as seen later. 

Before we study the region < A < 1 , we briefly look at the region — 1 < A < 
considered in ref. [17]. All the rapidities are symmetrically arranged on the real axis 
at $ = 0. We denote them by — oo < Ai < A2 < • • • < Am < 00. The largest 
rapidity Am goes to infinity at $ = 2(tt — [/,). Here the Am jumps from the real 
axis where a negative-energy mode lies to the in line where a positive-energy mode 
lies. This means a particle creation. When <3> increases from 2(ir — //), Am moves 
backward along the in line and the other rapidities remain on the real axis. The 
Am on in line may be called one-string. When $ = 2n, the rapidity Am reaches in. 
When $ = 27r — 2(7r — yu), Am goes to — 00 + in and jumps onto the real axis. This 
means a particle annihilation. The state comes back to the original ground state 
at $ = An. The relation between the initial and the final momentum and the one 
between the initial and the final phase shift are the same as in the regime A < — 1, 
which we will show in §4. The Berry phase of the ground state wave function in 
this regime has been calculated in ref. [18] to be n. 

In the case of half-filling M = N/2 which we consider, the picture of the Dirac 
sea emerges as follows. The eq. (1.5) reduces to pj = (Ij/N) x 2n with the 
condition $ = 0. The single-particle energy is Sj = — 2cospj = —2cos(2nIj/N). 
The negative energy modes are in the region < N/A and are given by the set 



M-l M-3 M-l , , 

h = jj— , ( 2 - 5 ) 



The same set {Ij} turns out to provide the ground state of the half-filling sector 
in the interacting case. 
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Figure 2-1. The momentum function and the phase-shift function for 
lm\ = (normal lines) and for Im\ = tc (dotted lines). 



In the next two subsections, we consider the case in which only a two-string 
is formed (0 < A < cos(7r/3)) and the case in which a three-string is formed 
(cos(7r/3) < A < cos(7r/4)). The conclusion we will draw in the former case is 
mainly based on the analytic study we have made at a finite size system M = 
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2, N = 4, which we present in Appendix A. The conclusion we will draw in the 
latter case is originally derived from the numerical study we have carried out in 
the case M = 3, N = 6. In both cases, we check the consistency condition that the 
momenta and the phase shift functions for the final configuration follow from those 
for the initial one by using the relations we derive from the global behavior of the 
rapidities. We find that this consistency condition is a nontrivial check to make. 
In the last subsection, we are able to see that the global behavior of the rapidities 
we find in the case of n-string formation ( cos(7r/n) < A < cos(7r/n + 1) ) passes 
this consistency check for any n. This check is very stringent and consolidates our 
conclusion. 

It will be helpful to point out here that one can predict a few things about the 
behavior of the rapidities without looking at the details of eq. (1.5). As we vary 
the twist angle continuously, an individual momentum, which is a root of eq. (1.5), 
must vary continuously. The only place in which a set of rapidities can jump to 
form a string while keeping the momenta continuous is the point at infinity ( i.e. 
the real part of the rapidity is infinite. See eq. (2.1)). This is also the place for the 
formation of a zero-energy bound state. ( See eq. (2.4).) The value of the twist 
angle where the formation of n-string at infinity can occur can be obtained from 
eqs. (1.5), (2.1). We start from eq. (1.5)for M rapidities and send the largest ( 
or smallest) n rapidities to plus ( or minus ) infinity, which subsequently form an 
n-string. Only when the twist angle $ = +(or— )n(n — (/), the rest of the M — n 
rapidities are kept symmetric with respect to the origin, staying in the ground state 
configuration. In other words, only at this value the ground state configuration for 
M rapidities reduces to that for the M — n rapidities. 

2.1. Ground state in the case of two string 

We now study the more specific case < A < cos(7r/3) (^ < /i < ^7r) where 
the two-string is expected to be formed during the process. For convenience we 
describe the process to be associated with the change of $ from — 2n to 2n. We 
first examine the latter half < $ < 2n of the process. When <3> increases from 
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zero the largest rapidity Xm reaches oo at $ = 2(n — fi) just as in the — 1 < A < 
case. The rapidity Xm jumps to the in line and begins to move backward along the 
in line from oo + in . The dramatic change arises when <£> increases further. The 
rapidity Am on the in line returns to oo + in. The second largest rapidity Xm-i 
on the real axis goes to infinity. The two rapidities Xm and Xm-i attract each 
other. At $ = 4(n — /i), the two rapidities reach Xm = oo + in and Am-i = oo. 
The two largest rapidities then jump together by in/ 2 to the imaginary direction 

O 1 

and they form a two-string: Xm = oo + ^in and Xm-i = oo + ^in. This means 
a formation of the bound state over the vacuum. The two rapidities are separated 
by in and their center is on the in line. So the deviation due to the finite size in 

(2) (2) 

(1.7) is given by 5m = S\ = ifi and 5m~i = ^ — 

When we further increase <3>, the two-string moves backward from plus infinity. 
The center moves on the in line from oo + in. In this process the difference 
between Xm and Xm-i decreases from in. At $ = 2n, the two-string arrives at 
the imaginary axis: 

X M = i(2n - /i), 

(2.6) 

Xm-i = 

We note that the relative rapidity {Xm — Xm-i) is 2i{n — fi), and the deviations 

(2) (2) 

5\ and 5^ in (1.7) vanishes. It is interesting that the momenta p(Xm), p(Xm-i) 
and the phase shift 9(Xm, Xm-i) are divergent at this point. The points (2.6) are 
the branch points of both the momentum p(A) and the phase shift 9(X), which are 
logarithmic functions (2.1) and (2.2) respectively. If the rapidity Aj goes around 
the branch point A = ifi (A = i(2n — fi)) counterclockwise, the momentum p(Xj) 
increases (decreases) by 2n. If the relative rapidity A goes around the branch point 
A = 2i(n — fi) counterclockwise, the phase shift increases by 2n. Therefore the 
behavior of the rapidities around these branch points has a global meaning and 
is necessary to determine the final configuration of the system. The paths of the 
rapidities Xm, Xm-i and those of the relative rapidity Xm — Xm-i around the 
branch points are plotted in Fig. 2-2 and Fig. 2-3 respectively. The Xm and Xm-i 
collide into the singularities perpendicularly to the imaginary axis (Fig. 2-2). The 
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value of the relative rapidity Xm — Am-i reduces from in to 2i(n — fi) (Fig. 2-3). 
This behavior of the rapidities is checked both by the numerical calculation and by 
the analytic calculation in the case of M = 2 which is described in Appendix A. The 
first half process in — 2n < $ < is similar to the latter half process. At $ = — 2n, 
the two-string starts from the singular points Ai = i(2n — n),\2 = if-i- When 
$ = — 4(7r — n) the two-string reaches Ai = — oo + i^n and A2 = —00 + i^. After 
jumping by —in/2, Ai and A2 move on the in line and the real axis respectively. 
Finally Ai jumps onto the real axis at $ = — 2(n — //). 



Im X 
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<2(7t-|i) 
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(a) the path of Xi 



Im X 




(b) the path of Xm 
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(c) the path of X2 



(d) the path of X M-1 



Figure 2-2. The paths of the rapidities in case of n/2 < ji < 2n/3. 



The relations between the initial state ($ = — 2n) and the final state ($ = 2n) 
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Figure 2-3. The paths of the relative rapidities in the case 7r/2 < ji < 2n/3. 



are determined from this behavior of the rapidities. They are rearranged as 

\(/) _ x (0 x (/) _ x(0 
(/) _ (4) (, = !,..., M-2). (2.7) 

Let us determine the relation between {p^, #(p^\pj^)} and {pP ,0{p[f\p^)}- 
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The cut of the function p(X) extends from the branch point ifi to another branch 
point i(2n — fi). The path of Ai starts from the branch point i(2n — fi) and that 
of Xm reaches the same branch point i(2n — fi) (Fig. 2-2). To determine the 
relation between p M ^ and p^\ we follow a circuit consisting of the path of Am, 
that of Ai and the real axis. We start from the branch point A^ = i(2n — fi) 
perpendicularly to the imaginary axis. Finally we come to the same branch point 
= ^ n ~~ Z 1 ) P er P en dicularly to the imaginary axis after running through the 
circuit. The circuit encircles the 'half of the branch point A = i(2n — /i) of the 
function p(X). We conclude that p(X M ^) is larger than p(X^) by tx: 

Pm = pf + (2-8) 

Similarly we find the relation 

Pm-i = P2 + t- (2-9) 

The other relations are trivial 

P { P=P% (j = l,...,M-2). (2.10) 

The cut of the function 0(A) extends from the branch point 2i(7r — fi) to another 
branch point 2%fi. To determine the relation between 9(p$ ,Pm-i) an d Q(P\\p$) 
we consider the circuit consisting of the paths of the relative rapidities Xm — Xm-i, 
Ai — A2 and the real axis. This circuit crosses this cut (Fig. 2-3). Therefore the 
relation between the two phase shifts is 

0(p$,P$-i)=8(p?,$)+to- (2.11) 

To consider the other relations for the phase shifts, recall that the difference in the 
two rapidities in the two-string is larger than 2(ir — fi) in this process. From this we 
say the relation Im (Am) > 71 + ij 1 — A*) = 2ix — fi and Im (A m-i) < 71 — i 11 — /-0 — f 1 
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in the two-string. Therefore Im (Am — Xj) > 2ir — fi > 2fi (j ^ M, M — 1) and 
Im (Xm-i — Xj) < /j, < 2(n — n) (j 7^ M, M — 1) are satisfied in the region of our 
consideration |r < /i < ^it. This means that the circuit consisting of the paths of 
Xm ~ Xj, that of Ai — Xj + 2 and the real axis does not cross the cut of the phase 
shift function (Fig. 2-3). Similarly the circuit consisting of the paths of \m-i~ Xj, 
A2 — Xj + 2 and the real axis does not cross the cut (Fig. 2-3). From these results 
we obtain the relations 

^ ) 4 ) -> = ^ ) 4U (212) 

where j = 1, . . . , M — 2. The other relations are trivial 



0(pfjP) = 0(p# 2 >p2.2)> ( 2 - 13 ) 



where j, k = l,...,M-2. 

The set {pf^} and the set {p^} are not equivalent even up to 2ir, as seen in 
(2.8) and (2.9). In spite of this, the initial state (<£> = — 2n) and the final state 
($ = 2n) are found to be the same. We show this in Appendix B. 

Our conclusion (2.8), (2.9) and (2.10) are different from the one derived in [19]. 
In ref. [19], it was concluded that the two-string never comes close enough to the 
iir line. They have obtained the relations p$ = Pi \ Pm-i = P2 an< ^ P^P = Pj+2- 
Their conclusion is not valid however. By summing over the relations (1.5) for all 
j (j — 1, . . . , M), we obtain a simple relation 

M ^ 

3=1 

which must be satisfied. Their conclusion does not satisfy this simple consistency 
condition. 



16 



Finally, let us check briefly the consistency between the boundary conditions 
(1.5)and our relations (2.8)-(2.13). The boundary conditions for the final state are 

2M P f + £ e(p { P,Pk ] ) = 2?r ( 2j ~f r ~ 1 ) + 2tt, (j = 1, . . . , M). (2.15) 
These equations are rewritten in terms of the initial set by using (2.8)-(2.13). 

2M(p\ 1 + n)+ ^fl(pf ) ,p? ) ) + 27rj = 2vr (^7T^) +2^, (2.16) 

2M( P « (X>(P« P?) - 2tt j = 2vr (^*) + 2vr, (2.17) 



(j = l,...,M-2). 
They are equivalent to the boundary conditions for the initial state 



(2.18) 



2M P % + yj <„rf ) ) = ^( 2fa " +2) 2 ~ Jtf ~ 1 )-^. 

V ^ / (2.19) 

(j = -l,0,...,M-2). 



2.2. Ground state in the case of three string 

Let us now turn to the case of three-string which is formed in the region 
cos(vt/3) < A < cos(vt/4), (g7r < jj, < j7i). We first describe the behavior 
of the rapidities we have found in this region. When $ increases from zero to 
2(vt — //), the largest rapidity Am goes to infinity and jumps onto the in line. 
When $ increases from 2(n — fx) to 4(vt — /i), the second largest rapidity Xm-i 
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goes to infinity and the two rapidities jump by in/ 2 simultaneously. Thus the two 

O 1 

rapidities Am = 00+ ^in and Am-i = 00+ ^in form a two-string as seen in §2.1. As 
$ increases further from A{n — /i), the two-string starts moving from plus infinity. 
The center of the two-string is on the in line and the distance between Am and 
Xm-i becomes shorter. The two-string and the third largest rapidity Am-2 attract 
each other. When $ = 6(n — //), the two-string and Am-2 go to infinity. Here 

the distance between Am and Am-i becomes 2n/3. Three rapidities are, therefore, 

4 2 
Am = 00 + g27r, Am-i = 00 + g27r, Am-2 = °0- Here the three rapidities jump 

by ^ji. They form a three-string where Am = 00 + ^in, Am-i — 00 + in and 

Am-2 = 00 + \in. The deviation of the length of the string due to the finite size in 

(1.7)is 5 M = S[ 3) = 2i(§7r - //), 5 M -\ = 5^ = and 5 M ~2 = ^ = — 2i(§7r - fi). 

When $ increases from 6(n — /i) to 4n — 6(n — fi), the center of the three-string 
moves from 00 + in to —00 + in. The three-string fluctuates in this process (Fig. 
2-4). It is interesting to note that the deviation 5m(^m-2) moves counterclockwise 
(clockwise) M — 1 times around the origin (5 = 0) which is the singular point where 
the phase shift diverges. This singular point in the relative rapidity space is the 
branch point of the phase shift function #(A). Thus how many times the relative 
rapidity winds around the branch point determines the corresponding increment 
of the phase shift. When the relative rapidity A moves counterclockwise around 
the branch point A = 2i(n — /i), the phase shift increases by 2n. The phase 
shift 9(Xm, Am-i) (#(Am-i, Am-2)) increases (decreases) by 2(M — l)n due to the 
fluctuation of the three string. 

At $ = An — Q(n — /i), the center of the three-string reaches —00 + in. The 
rapidities jump by — ^ and are decomposed into Am-2 on the real axis and the 
two-string (Am — Am-i) whose center is on in line. At $ = An — A(n — /i) the 
two-string is decomposed; Am and Am-i jump onto the in line and the real axis 
respectively. Finally Am jumps on to the real axis at $ = An — 2(n — yu). When 
$ = An, the set of the rapidities {AV '} is the same as the set of initial ($ = 0) 
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Figure 2-4. The fluctuations of the three-string in the case of M=3 and 
N=6. 



rapidities {A^}. The rapidities are rearranged as 

A (/) _ x (0 A (/) _ x (0 x (f) ... x (0 
Af = A«3(j = !,•••, M-3). 

This behavior of rapidities has been checked by numerical calculation (Fig. 2- 
5). The motion of the two-string is independent of the rapidities lying on the real 
axis as seen in § 2-1. In the case of the three-string, its fluctuation is influenced 
by the rapidities lying on the real axis. 

The relations between the set {pf\ 9(p^/\pi^)} and the set {pf , Q(j^j\Pk)} 
are determined from the behavior of the rapidities. The branch cut of the function 
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. 3- 2n I 3 - 2(7t-|i) 



0.25 - 
0.2 - 

. 15-- 
. 1- 

. 05- 



Figure 2-5. The deviation £3 = (A3 — A2) — 2(n — li) in the case of 
M = 3, N = 6 and /i = 5n/7. When $ = 6(n — li), the deviation is 
2n/3 — 2(7r — li) ~ 0.299. When $ increases from G(ir — fi) to 47r — 6(7r — fx), 
the deviation moves counterclockwise twice around the origin. 

p(X) extends from the branch point A = ifi to another branch point A = i(2ir — fi). 
Only the path of A^-i crosses the cut (Fig. 2-6), so that the relations of momenta 
are 



The cut of 6>(A) extends from the branch point A = 2«(7r— /i) to another branch point 
A = 2i(i. The path of the relative rapidity \m — Am-1j an d those of \m-i ~ Am-2 
and Ajvf — \m~2 cross the cut 1 + (M — 1) times, 1 — (M — 1) times and once 
respectively (Fig. 2-7). We have, therefore, 



-0.1 -0v0.5 

-0 . 05- 



. 1 



Re5 3 




(2.21) 




(2.22) 
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Im A, 




(b) the path of A, M -i (c) the path of I 



Figure 2-6. The paths of the rapidities in the case of 2n/3 < /i < 37r/4. 

For the other phase shifts, only the path \m-i — Afc (A; = 1, . . . , M — 3) crosses the 
cut once (Fig. 2-8). We have 

^( A M > 4 /} ) = A fc+3)> 
^ ) _ 1 ,Ai /) ) = ^ ) ,A« 3 ) + 2vr, 

6(\$_ 2 AP) = e(>$\\% i ), (2-23) 
^(Af,A^) = ^3,A« 3 ), 

(j,A; = l,...,M-3). 

Finally, let us check the consistency between the relations (2.21)-(2.23) and 
the twisted boundary conditions (1.5). The boundary conditions for the final state 
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(b) the path of (k MA - X M _ 2 ) 



(c) the path of (^ M -^ M _ 2 ) 



Figure 2-7. The paths of the relative rapidities in the case of 2n/3 < /x < 
3tt/4. 



are 



2j - M - 1 



+ 4tt, 



(2.24) 



(j = l,...,M). 

These equations are rewritten in terms of the initial set by using (2.21)-(2.23), 



M — 1 



4vr, (2.25) 



2M(p^+2tt)+ [ 9 (P2 iPjU) - 2M7r - 2 ( M - 2)vr + 2(M - 3)tt 



2tt| ^-JL ]+4vr,| 
(2.26) 
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(a) the path of (X M -Xj) 

0=l,2,...,M-3) 



Im A, 
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2(7t-n)> 




' 





Im X 



2n 

.2(tc-h) ; 



Re X 



(b) the path of (V-a-i) 

(j=l,2,...,M-3) 



(c) the path of (X M . 2 - Xj ) 

0=l,2,...,M-3) 



Re A 



Figure 2-8. The paths of the relative rapidities in the case of 2n/3 < \x < 
3tt/4. 



2Mpf + [ Y. 9 ^i\pf+z) + 2 ( M - 2)tt - 2vr ] = 2tt 

,fc^3 



M-5 



+ 4tt, (2.27) 



2- [ 2j M 1 )+4tt, (j = 1,...,M-3). 



(2.28) 



23 



They are equivalent to the boundary conditions for the initial state: 



2M Pl + J20(P i i\p ( i ) ) 



k=£3 



Jfc^l 



k=£2 




+ 0, 



+ o, 



+ o, 



(2.30) 



(2.31) 



(2.29) 



2Mp j+3 + J2 e(p%i,p®)=2K 




) 



+ o, 



(2.32) 



2.3. Generalization to the case of ^-string 

We have found the collision of the rapidities forming the two-string into the 
singularities as well as the fluctuation of the three-string around its center in the 
previous subsections. It is possible to generalize these behaviors to the case of 
n-string. Let us describe the motion of the rapidities we have found in the region 
cos(^) < A < cos(^py) {^n < n < ^n). At $ = 0, all rapidities arc 
symmetrically arranged on the real axis with respect to the origin. At <£> = 2j(n—fi) 
for (j = 1,2, ■ ■ • ,n), a group of j rapidities at Re(X) = oo jumps: a j-string is 
formed from a (j — l)-string and one string . The outcome of this series of events 
is the formation of the n-string at $ = 2n(n — u). After this, the center of the 
n-string moves backwards on the in line from oo + in to — oo + in as $ increases 
from 2n(7r — fi) to An — 2(n — fi). The n-string itself fluctuates during this period 
as we have seen in the case of the three-string. (We will describe this shortly for 
the case n even and for the case n odd separately.) As $ increases further, a 
separation of a j-string into a (j-l)-string and one string on the real axis occurs at 
$ = 47r — 2j(n — /j,) for j — n, n — 1, • • • , 1. Finally at $ = An, all rapidities are 
again symmetrically arranged on the real axis with respect to the origin. 
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Let us now describe the fluctuation of the n-string. In the case of n even, the 
two rapidities in the middle, namely, \m~^+i and \m~% of the n-string collide 
into the two branch points of p(X) at $ = 2n: 

Ajvf-f +1 = 27r - u, \m-% = (J- (2.33) 
The other rapidities fluctuate: the deviation 

Tl 

5 M - j+ i = (Xm-j+i - A M -j) - 2(tt - ft) (j = 1, 2, 1) (2.34) 

moves counterclockwise (M — n + 2j) times around 5 = 0, which is the branch 
point of 0(A). The deviation 

$M-n+j = (^M-n+j ~ ^M-n+j+l) ~ 2(7T - //) (j = 1, 2, 1) (2.35) 

moves clockwise (M — n + 2j) times around 5 = 0. 

In the case of n odd, the rapidity in the middle A M _n^i moves on the m line 
while the other rapidities fluctuate: the deviation 

Tl — 1 

f>M-j+i = (Ajif-j+i - Xm-j) - 2(7r - ix) (j = 1, 2, • - • , — ^ — ) (2.36) 
moves counterclockwise (M — n + 2j) times around 5 = and the deviation 

Tl — 1 

$M-n+j = {^M-n+j ~ ^M-n+j+l) ~ 2 (7T - fj) (j = 1,2, - ) (2.37) 

moves clockwise (M — n + 2j) times around 5 = 0. Only the rapidity A Af n-i 
crosses the cut. 

The rationale for this behavior of the rapidities lies in the fact that they pass a 
set of consistency conditions, which is found to be very stringent. Let us denote by 
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{A^} a set of initial rapidities at $ = — 2-k and by {A^} a set of final rapidities 
at $ = 2n. We have 



A 



(/) -A« (, = 1,2, 



M-j+l ■ -j 

x (/) _ x « 
A 3 ~ A j+n 



(j = l,2,---,M-n) 



(2.38) 



From the paths of the rapidities described above, we obtain a set of relations for 
the momenta 



.(/) 



(0 



ra — 1 



p(A« i+I )=p(A<")+-, P(A^U)=P(AP +1 )+- (-even). 



P( X M-^i) = P ( A 2±) + 2tt (n : odd). 

1V1 2 2 



p(A^)=p(A^) (j = l,2,---,M-n). 



(2.39) 



From the paths of the relative rapidities, we obtain a set of relations for the phase 
shifts: 

8&M-i + v *M-i) = d ( X f> >&i) + 2 ^{1 + (M - n + 2j)}, 



(.7 = 1,2, 



n — 1 



(2.40) 



^i?-^* = » (Ajjl i+1 , Af ) + 2vr{l + (M — n + 2j)}, 



(J = l,2, 



n — 1 



(2.41) 



0(a{£_» +1 , ^m-s) - ^(An^p Al v ) + 2tt, 



>(/) 



(0 



X M-k+i) = AH + 2^ (1 < J < h < n and * - j > 2), 



.(0 \« 
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*(*m > Al /} ) = 9(X?, A«J (fc = 1, • • • , M - „), 
'(4?-,H-i» 4 /} ) = ^J?, A« J (k = 1, • • • , M - n), 
Q^M-j+v Ai /} ) = 0(*f , Ag. n ) + 2vr (2 < j < n - 1, fc = 1, • • • , M - n), 

Oixf, \ { k f) ) = 9(\fl n , A« J (j = 1, • • • , M - n, k = 1, • • • , M - n). (2.42) 

In the first and the second formulas, the factor ±(M — n + 2j) represents a contri- 
bution due to the fluctuation of the n-string. The third formula is for the case of 
n even only. Following the same procedure as described in §2.1 and§2.2, and using 
the relations given above, we have checked that the twisted boundary conditions 

(i) 

for the set of the initial rapidities {A^ } follow from those for the set of the final 
rapidities {A^ }. The proof is too long and tedious to be presented here. 

This consistency check is in fact very stringent. We are convinced that the 
picture of the n-string fluctuations we have found in this paper has now been 
established by this. 
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3. Results in the Regime A < —1 

In this section, we study the ground state properties mentioned in the intro- 
duction of the XXZ model in the regime A < — 1. Throughout this regime, the 
single particle energy is kept non-positive. 

3.1. THE REGIME A < -1 

In this regime, the momentum is parametrized in terms of the rapidity as 



p(X) = — iln 



sinh^(*A + /i) 
shrill (iA — /i) 



(3.1) 



where A = — cosh//. As the momentum increases from zero to it (from n to 2n), 
the rapidity increases from zero to n (from n to 2ir). The particle-particle phase 
shift becomes 

sinhi(i(Ai - A 2 ) + 2/x) 



0(Ai,A 2 ) =iln 



sinh^(i(Ai — A 2 ) — 2/x) 



(3.2) 



As the relative rapidity increases from zero to n (from n to 2n), the phase shift 
decreases from zero to —71 (from — 7r to — 2n). We plot the momentum p(X) as 
a function of the rapidity and the phase shift 0(A) as a function of the relative 
rapidity A in Fig. 3-1. The total energy is expressed as 



A 7 



E = — cosh/i + ^ 



— 2sinh [i 
J coshyU — cosAj 



(3.3) 



The single-particle energy 



£ j = 



-2sinh ji 



cosh/i — cosAj 



(3.4) 



is negative for the real Xj. 



28 




Figure 3-1. The function for the momentum and the function for the 
phase-shift function in the case of A < — 1. 

Let us consider the process associated with the change of the twist angle <£> from 

zero to 47r. The motion of the rapidities is determined from equation (1.5) with 

(2.5). When $ = 0, all the rapidities Xj are real and symmetric with respect to the 

origin of the rapidity space {{Xj} = {— Xj}) within the Brillouin zone — n < X < it. 

At $ = 27r, the largest rapidity Xm (or one-string) reaches 7r and the rest of them 

are symmetric with respect to the origin of the rapidity plane A = 0. This is an 

excited state for the periodic boundary condition. At $ = 4n we will denote the 

( f) 

configuration of Xj by {A^ }. It is related to the original ground state configuration 
by 

A<£> = A< 4) + 2tt, Af = A« 1 (j ^ M). (3.5) 
This leads to the relations 

pff = pl° + 2tt, =pgi U*M), (3-6) 



29 



0(p i I ! I \p ( j f) ) = O(p?,pf +1 )-2n, (3.7) 

0(pfA f) ) = °(pf + vPkli) U ^M,k^ M), 

where p^ = p(X^) and = p(X^). The set of the momenta for the final 
state coincides with the one for the initial state up to 2n. Accordingly the system 
returns to the ground state at $ = An. 

We check briefly the consistency between the above equation and the twisted 
boundary condition (1.5) at $ = 0,27r, 4n. The twisted boundary condition for 
$ = is 

2Mp{xf) + e(\f, xf) = 2n 2j -f~\ (3.8) 

k+j 

As the p(X) and 0(A) are odd functions, the above equation becomes 

2M P (-Xf ) + £ 6(-Xf\ -A») = 27r -M + 2(Af-i + l)-l j (3g) 

Comparing this with the twisted boundary condition for p{X l M _- +l ) 
2Afp(AS_ J+1 ) + ^ g(A^ _,- +1> A* M _ fc+1 ) = 2tt ~ M + 2(M ~ J + 1} ~ 1 , (3.10) 

(i) (i) 

we find that the configuration of the rapidities is symmetric Ay_ J+1 = — A!- . In 

the case of $ = 2n, we can derive A^ = n and Xm-j = —Xj (j ^ M) in the same 

f f) 

way. When $ = 47r, the twisted boundary condition for p{X - ) is 

2Mp{xf) + ^ g(Af , A^) = 2n 2j " ^ " 1 + 4tt. (3.11) 

Using Eqs. (3.6) and (3.7), the left hand side of the equation is rewritten in terms 
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ofAf, 

2M P {\f +1 ) + ( 0(*? + v A?) + A = 2* 2j -™- 1 + 4n, (3.12) 
for j ^ M and 

2M P (A«) + 4°) - 2-) = 2rc^ + 4vr, (3.13) 

where we used the relation 9(Xj,Xk) = — 9(Xk, Aj). Eqs. (3.12) and (3.13) are 
equivalent to the boundary conditions for the initial state 

™ P (X%) + 2 Xf) = 2, ^ + 1 »- M - 1 + 0, (3.14) 

2Mp(A«) + ^ 6>(A« A<<>) = 2. (-^) + 0. (3.15) 

We confirm the motion of the rapidities for the case of M = 2 and iV = 4 in 
Appendix A. Observe that the single-particle energy (3.4) remains negative during 
the process. 

3.2. THE A = -1 CASE 

Let us now study the ground state in the critical coupling case A = — 1. When 
A = — 1, the momentum p is parametrized in terms of rapidity A by 

p(X) = 2arctan(2A). (3.16) 

When the momentum increases from zero to n, the rapidity goes from zero to +00. 
As the momentum increases still more the rapidity jumps from +00 to —00. When 
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the momentum increase from it to 2tt, the rapidity goes from — oo to zero. The 
phase shift and total energy become 

0(Ai,A 2 ) = -2arctan(Ai - A 2 ) (3.17) 

and 

N M -1 

E = - + Y (3.18) 

j=i J 

respectively. The single-particle energy 

-1 

£j ~T74a] 

is negative for real Xj and zero for Xj = ±oo. 

Let us describe the process associated with the change of twist angle <3>. When 
$ = 0, all the rapidities Xj are real and symmetric with respect to the origin of 
the rapidity space. If we set $ = 2n, the largest rapidity Xm goes to plus infinity 
and jumps to minus infinity. At $ = 4tt the configuration of rapidities {Xj} comes 
back to the initial configuration as follows 

Ai? = A« A^ = A« 1 (j^M). (3.20) 

This leads to the relations 



(3.19) 



p { £> =pf + 2tt, pf=p% (j^M), (3.21) 



e(pZ\p ( ? ) ) = e(p?4l 1 )-2n, 

This is the same as the relations (3.6) and (3.7) in the A < — 1 case. The period 
of the system is again A$ = An. 
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We have found that the motion of the rapidities, the period, the final config- 
uration of the momenta and that of the phase shift functions in the case A < — 1 
are the same as those found in the case — 1 < A < 0. We have, therefore, es- 
tablished the continuity of the two regimes. This can be heuristically explained 
by the structure of the degeneracies of the levels, which we illustrate for the case 
M = 2, N = 4 in Fig. 3-2. The ground state at $ = never gets degenerated with 
the zero energy levels as long as A < 0. This is, however, seen for A > 0. 




(a) A<0 (b) A=0 (c) A>0 



Figure 3-2. The energy levels of XXZ model in M=2 and N=4 case. The 

eigenenergies are E = A ± y/A 2 + 8 cos 2 $/4, A ± ^A 2 + 8sin 2 5/4, 0, 
which we derive in Appendix A. 
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4. The Berry Phase 



Let us now turn to the discussion of the Berry phase associated with the 
twist angle $(£), which is now considered to vary with time. We are interested 
in computing the nontrivial geometrical phase 7 appearing in the wave function 

M*(t),f)> = exp j-i ^oW))^ j exp(i 7 (i))|VoW))>- (4.1) 



The substitution of this into the Schrodinger equation provides a formula 

t 



7(t) = Re 



i /(it 



(^o($(t))|^o($(t))) ^ 
(V>o(S(t))|^(S(t))) 



(4.2) 



In our consideration, the initial and the final twist angles are taken to be $(0) = 
— 2n and <£>(£) = +2n respectively. After a change of variables, we find 

2tt 



7 = Re 



L -2tt 



(M*)\&\M*)) 
(M*)\M*)) 



(4.3) 



We compute this quantity in different regions. 



4.1. A < -1 



We take the half-integer coordinate Xj = — ^f^, ■ ■ ■ , where N = 

2M. We first examine the relation between the wave function x(^0 °f the initial 
state at $ = — 2n and that of the final state at $ = 2ir. We first want to establish 
( see (1.2)) 

| X (-2tt)> = -|x(2tt)>. (4.4) 



N-l 



By using the permutation 



Q'a = 



Q a + 1 (for Qa^M), 



(4.5) 



1 (for Q a = M), 

and the relations (3.5), (3.6) and (3.7) we find the difference in the phase factors 
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from the initial state to the final state 



(-l)Wi = (-l)ivi x (-1) 



\Q'\ 



xM-l 



(4.6) 



M \ / M 

exp [iJ2 x *P { Q a = ex P pX^gl ) X 



(4.7) 



a=l 



a=l 



ex p ( ££*(p$>p£m*6-*«) ) = ex p ( IE^S'pS)^-^) ) x 

a<b 



M-l 



a<b 



(4.? 



Thus under the full period $ = 47r, x gets multiplied by a factor (— 1) M x (— 1)(— 1) M 1 
— 1, which is what we wanted to show. 

Identifying 

|^($))=exp(-i$/4)|x($)), (4.9) 

We find 

2tt 



7 = 7r + Re 



V X<yXh XM )*fe( Xl ' 111! 
■ x(^b • • • , ^m)*x(^i, • • • , xm) 



. -2?V i X i} 



(4.10) 



The summation represents fl^i where each runs over the region of 

half integer values. The function of <3> 



{Xi} 



X(xi, ■ ■ ■ , x M \pi($), ■ ■ -Pm(^))*^x(xi, x M \pi($), • •■Pm(^)) 

IIxWII 2 



(4.11) 



is odd as shown in Appendix B. The second term in (4.10) is zero and the Berry 
phase is equal to 



7 = 7T. 



(4.12) 



Recalling the calculation of ref [18], we find that the Berry phase is n through- 
out the region A < 0. This confirms the continuity at A = — 1. It will be different 
in the region A > 0, which we will find in the next subsections. 
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4.2. 0<A<costt/3 



The strategy for the computation is the same. By using the permutation 



Q a + 2, (forQ ^M-l,M), 
Q' a =\ 2, (forQ = M-l), 

1, (for Q a = M), 



the relations (2.8)-(2.13) are expressed as 



(4.13) 



(4.14) 



e (Po!>PQ.) = 9 (PqvPqO + ^ 5 Q'^ 5 Q' a ,2 ' 2^ Q ,, 2 5 QU . (4.15) 
Counting the change of the phase factors 

(-l)IQI = (-1)10'! x (-1), (4.16) 



exp 



' M 
a=l 



,(/) 



exp 



M 



EM 



a=l 



exp 



a<b 



exp 



^0(Py b ,Pcy a Wx b -x a ) 

a<b 



(4.17) 

x (-1), 
(4.18) 



we find, in the same way as we have derived Eqs. (4.4) , 



|x(-2tt)> = |x(2tt)>. 



We conclude 



7 = 0. 
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4.3. cos 7r/3 < A < cos 7r/4 



The corresponding formulas in this case are 



Qa 



' Qs + 3, 


(for Q a 


^ M - 


2,M-1,M), 


3, 


(for Qa 


— M — 


2), 


< 


(for Qa 


— M — 


1), 


1, 


(for Qa 


= M), 





(4.19) 



(4.20) 



and 



+ 2(M - 2)7t(5q/ j2 %,3 - S Q ' b ,^Q'a,2) 
+ M 5 Q' b ,l 5 Q' a ,3 - $Q> b ,3$Q' a ,l) 

{MM ~\ 
S Q' b ,2 Yl 5 Q'»,k - Yl 6 Q' b ^Q' a ,2 \ , 



(-1)101 = (-l)IQ'l x (_l)(A / /-l)+(M-2)+(M-3) ? 



exp 



Elf) 



a=l 



cxp 



M 



(0 



a=l 



x (-1), 



(4.21) 



(4.22) 
(4.23) 



e 5Ea<»»(Pg>oi) £ (^-^) =e 3E.<6«(p2[j'2) l W^-^) x ( _ 1) M+(M-2) + l+(M-3) > 



We conclude 



7 = 7T. 



(4.24) 
(4.25) 
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4.4. cos(7r/n) < A < cos(7r/n + 1) 

By using the relations we have found in §2.3., we compute the Bery phase in 
this regime in the same way as we did in n = 2, n = 3. The final answer is 

7 = for n even 
7 = 7r for n odd. 

4.5. Discussion 

Finally we would like to argue that the Berry phase we have computed can 
be used as a measure of the statistics of the quasiparticle ( or the bound state) 
involved in the process. The above calculation shows 

7/vr = m mod 2. (4.26) 

for the case in which the m-string goes around the loop in the momentum/rapidity 
plane. This quantity 7/71 mod 2 measures the monodromy property in the p space 
of the particle involved. We propose to use this as a definition of statistics. To 
make our discussion a little more concrete, let us consider the case A = — 1. At 
$ = 27r, the largest rapidity ( or one-string) goes to the edge of the Brillouin zone 
( or infinity in the A plane.) This root ( or the rapidity) decouples from the rest of 
the roots and the system of M equations eq. (1.5)reduces that of M — 1 equations. 
According to the description of [21], this is a two-spinon state over the vacuum. 
Eq. (4.26)provides 1 for a creation and a subsequent annihilation of two spinons, 
which implies 1/2 for a spinon. 
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6. Note Added 

After the completion of this paper, we became aware that the formation and 
the collision of the two-string roots as a function of the coupling and the twist 
angle is mentioned by Alcaraz, Barber and Batchelor in [22]. We thank Professor 
Batchelor for informing us of this reference. 
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APPENDIX A 

In this appendix we examine the case M = 2, N = 4. The boundary conditions 
'1.5) are 

4pi -9{p2,pi) = -7r + $, (A.l) 



4p2-0(p 2 ,Pl) = +7T + $, 



(A.2) 



while the phase shift 



9(p2,Pi) = 2arctan 



Asin[(p 2 -pi)/2] 



cos[(p 2 +Pi)/2] - Acos[(p 2 -Pi)/2] 



(A.3) 



In what follows, we solve explicitly this system of equations. The sum (A.l) + (A.2) 
leads to the relation 



P2+P1 



(A.4) 



So the phase shift (A.3) becomes 
0(p2,pi) = 2arctan 



Asin[(p 2 -pi)/2] 



cos[$/4] - Acos[(p 2 -pi)/2] 



(A.5) 



From difference (A.2) — (A.l) we obtain another formula for the phase-shift. 



0(P2,Pi) = 7T - 2(p 2 - pi). 



(A.6) 



These formulas of the phase shift (A.5) and (A.6) lead to the equation on the 
relative momentum 



2arctan 



Asing 



cos($/4) — Acosg 



ii — 4g 



where q = (j>2 — Pi)/2. This is expressed as a relation 



$ 2 $ 
2 cos —cos q — A cos o — cos — = 0. 
4 4 



(A.7) 
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When cos($/4) = 0, the equation becomes 



cos q = 0. 



When cos($/4) ^ 0, the roots of the equation (A. 7) are 



cosg = 



A ± A 2 + 8cos 2 ($/4) 
4cos($/4) 



(A.8) 



In the case | cosg| < 1, p\ and p2 are real number. In the case of | cosg| > 1, p\ 
and p2 are complex number with pi* = P2- This set of the complex momenta is 
the two-string. When cos($/4) goes to zero, the relative momentum becomes 

(+ ) r A + v/A2 + 8cos 2 ($/4) 
lim coso v ' = lim — — 

cos($/4)-»+0 cos($/4)->0 4cos($/4) 

+ oo for A > 0, 
for A < 0, 



lim cosg*- 

cos($/4)->+0 



lim 

cos($/4)^0 



A - y/A 2 + 8cos 2 ($/4) 
4cos(<3>/4) 
for A > 0, 



+ oo for A < 0. 



From (A. 4) and (A.8), p\ and p2 are determined as 



$ 

pi — - — arccos 
4 

$ 

P2 = — + arccos 
4 



A± v / A 2 + 8cos 2 ($/4) 
4cos($/4) 

A± v/A 2 + 8cos 2 ($/4)' 
4cos($/4) 



The wave function is determined from the momenta as follows. In the present 
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system the wave function (1.2) with (1.3) is 

\x) = Yl Yl e ^ 2 ~~ Xl "> 5Z( _1 )' QI ex P( ix iPQi + ix 2PQ 2 

xi x 2 Q 

x exp(-6{pQ 2 ,p Ql )e{x2 - xi)) 

IT) 

The wave function is rewritten as 



|X> = Y i ex P( ix iPi + ix 2P2 + ^0(Pl,P2)) 

{x 1 ,X 2 \x 2 >X 1 } 

- exp(ixip 2 + ix 2 pi - ^0(pi,P2))}cr~ l a~ 2 \ |) 
= 4i Y exp(i Pl ^ P2 (x 2 +a?i)) 

{x!,X 2 \x2>Xi} 

x sin ^ 2 - Pl (x 2 -an) + ^(pi,P2))<r-cr- 2 | |> 
By using (A. 4) and (A. 6), the wave function becomes 
\ X ) = U exp(i^(x 2 + x 1 ))cos( P2 2 Pl (g 2 -ari-2))<7~(7~l t)- (A.9) 

{xi,a;2|a;2>a;i} 

We take the coordinate as 



11 1 2'2'2'2' i 



The wave function is described in terms of the relative momentum; 

\ X ) = 4i]exp(-i|) | - | ^> + exp(+i|) | - ^> 

+ icos^— ^{exp(-4)| - ~) + I - ^, 1) + exp(+4) l> + ~ 



We get the explicit formula of the wave function by substituting (A. 8) into this. 
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The eigenenergy is calculated from (A. 4) and (A. 8) 



E = At a/A 2 + 8cos 2 ($/4). 



The energy, relative momentum and the momenta of the ground state are 



E<® = A - y^A 2 + 8cos 2 ($/4) 



(0) A + v/A 2 + 8cos 2 ($/4) 
cosg w = 



4 cos 



(o) 
Pi 



(0) 
P2 



arccos 

2 



— + arccos 



A + y/A 2 + 8cos 2 ($/4) 



4 cos * 



A + v/A 2 + 8cos 2 ($/4) 



4 cos 



Let us consider the case A < 0. In this case the momentum is real number 
because of 

A + y/A 2 + 8cos 2 ($/4) 



cos q 



(0)1 



4 cos j 



< 1 



The motion of the momentum in the process is as follows. When $ = 0, the 
momenta are 



(0 

p\ ' = —arccos 



P2 — +arccos 



A + VA 2 + 8' 



A + VA 2 + 8' 



(A.10) 



thus the momentum is arranged symmetric pp = . When $ = 2ir, cosq^ = 
0, so the momenta are 

Pi = 0, P2 = 7T. 

When $ = 47r, the state comes back to the ground state. The set of momenta 
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{P'liPl} at $ = 47r is equal to the set of momenta {Pi'jP^'} at $ = up to 2tt; 

f 4tt A + VA 2 + 

p\ = — arccos 



,(0 JOi 



7r — 7r — arccos 



A + 



—arccos 



(») 

= P 2 , 



A + VA 2 + 8' 



f 47T 

p 2 = + arccos 



A + VA 2 + 8' 



= 7r + 7r — arccos 



A + VA 2 + 8' 



(A.H) 



Our assertion on the motion of the rapidities in Sec. 4-3 is checked by (A.lO)-(A.ll). 

Next we consider the case A > 0. We take the region — 2tt < $ < 2n for 
convenience. Note the relation 



cosg 



(0)1 



A+ v/A 2 + 8cos 2 ($/4) 



4 cos 



<1 for -2(n-fi) < $ < 2(n-fi), 
> 1 for - 2tt < $ < -2(7r - /j>), 2(n - fi) < $ < 2n, 

we check that the two-string is formed at $ = 2(n — fi). In the two-string the 
momenta is complex number with p\ = p 2 which is equivalent to the center of 
rapidities (Ai,A2) is on the in line (see (2.1)). When $ goes to 2tc, the relative 
momentum become — > too because of cos — > oo, so the momentum becomes 



p{ 
f 

P2 



71 
7T 

2 +l °°- 
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From this and Eq. (2.1), we find the rapidities arrived at the singular point 



Af = i(27r-A*), 
\{ = ifi. 

In the same way, we obtain the momenta and rapidities at $ = — 2n as 

i n 
Pi = + ^ 00 ' 



P2 



71 

ioo, 

2 



So, we get the relations 



\\ = ifi, 

\ % 2 = i{2ir - fj). 



\f — X* 
A 2 — A l5 



\f — X* 
A l — A 2i 



P2 = Pi + 7T, Pi = P\ + 7T- 

They agree with (2.7)-(2.8). 

The set {pj} is not equivalent to the set {p 1 -} up to 2tt. However the initial 
state and final state are the same state. The relative momentum cos = oo at 
$ = ±2n, thus the wave function (A. 9) becomes 



|x(±2?r)) oc |exp(-i 



±2tt, 



3 1, , ±2vr, 

2'2> +6Xp( ^ 



13 | 



+ i 2 oo x ( exp(— i 



±2tt, 



3 1, 11. , ±2vr 1 3 X 

-2>-2> + !-2'2> + eXp( ^ )l 2'2> 



3 

'2' 



the normalized wave function is 
Ix(±2tt)> : 



. .±27T.. 3 l x 

ex P H^-)| --,--> + 
o' o/ + I o' o) 



1 2' 2' 1 2'2 
Therefore we obtain |x(27r)) = |x(— 27r)). 



1 1. . ±2tt 1 3 x ,3 3 X 

-2'2> +eXp( ^ }l 2'2> + l 2'-2> 
1 3v 3 3," 

'2' 2' + l 2'~2/. • 
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Finally we examine the energy spectrum by the explicit diagonalization of the 
Hamiltonian (1.6) with the twisted boundary condition which is equivalent to the 
Hamiltonian 



(A.12) 

with boundary condition 

C'j+iv = Cj, Cj+iv = ^j' 

where 

The diagonalization of the Hamiltonian (A.12) for iV = 4, M = 2 system leads to 
the energy spectrum 



E = A± yA 2 + 8cos 2 (-), A± yA 2 + 8sin 2 (-), 0, 0. 

We plot the spectrum as a function of $ in the case A < 0, A = and A > in 
Sec. 4-4. 
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APPENDIX B 



The set {p^} and the set {p^} are not equivalent even up to 2n, as seen in 
(2.8) and (2.9). In spite of this, the initial state ($ = — 2ir) and the final state 
($ = 27r) are the same. To show this, we see that the wave function \x(~ 2tt) > 
of the final state (at $ = 2n) is invariant under the set of the transformations 

Pm + 71 and Pm-i Pm_i + n - When $ = 2n - e (0 < e < 1), the 

ff) (f) 
rapidities of two strings are \^ ~ i(27r — jj) + 5 and A)^-^ ~ i/i + 5, where 

< 5 <C 1. Substituting this into (2.1), we obtain 

P 2>~-iA + |, pg^-iA+l, (B.l) 

where A = ln(2sin///5) 3> 1; A goes to infinity for e — > 0. The phase shift 
%M'Pm-i) * s determined from the twisted boundary conditions (see (1.5)), 



2^^! + eipfl^pfl) + £ eipfl^pW) = 2rr {^j^j + 2vr. (B.3) 
These lead to the relation 

M-2 

2M(4 } -4U) + 26{pl\p { S_ 1 ) + ^(Pm^P) ~ Oipfl-vPP)) 

k=i 

=2ir. 

(B.4) 

From (B.l) and (B.4) the phase shift becomes 



e(P { £,P { £-i) ^ ^MA + vr - \ (0(p { £,P { k f) ) - 0(p$-vPk f) ))- ( B - 5 ) 

fc=i 

If we substitute this into the Bethe's wave function (1.2) and (1.3), the A depen- 
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dence of the wave function of the final state (at $ = 2n) becomes 
Ix(2tt)> 

= Yi n^-^a)!^- 1 ) 1016 ^ * 00 ^ 

xi,—,XMa<b Q j 

= E II e ( Xb ~ Xa ^ $^(- 1 )' Q ' ex P [ A ( x Qm - x Qm-i ~ M 4 x Qm ~ x Qm-i))] 
xi,—,XMa<b Q 

x (independent part of A) J^Jcr~ | f)- 

(B.6) 

Considering the range of the coordinate, it is shown that the leading terms (~ 
e A(M-i)^ Q £ ^ e wave function (B.6) satisfy the condition xq m — XQ M _ t = — 1 or 
2M-1: 

N-l _N-3 N-3 N-5 
\ x Qm-> x Qm-i) ~ \ 2 ' 2 2 ' 2 ' ' ' 

J-5 iV-3 JV-3 iV-1 JV-1 iV-1 ( B ' 7 ' ) 
" ' ' ' "I - ^' ' "I - ^' ' 2~~ 

These xq m and xq m1 are adjacent to each other on the ring. When e — > 
(A — > oo), the leading terms go to infinity faster than the other terms. Therefore, 
after the normalization of the wave function, the leading terms remain and the 
other terms vanish. The wave function at $ = 2n (e — > 0) is the sum of the 

terms which satisfy the condition (B.7). Under the set of the transformations 

C f ) ( f) ( f) ( f) 

Pm ~^ Pm + 71 an d Pm-i ~^ Pm-i + 71 ■> ^ e phase factor expi Yl x Q a Pa m the wave 

function is transformed as 

e iJ2 x Q a p* e *X>Q a Pa x e i7r ^ QM+:EQM - 1 ^ (B.8) 

As xq m + xq m1 is even as seen in (B.7), the phase factor in (B.8) becomes unity: 

Q i<xq M +XQ M _ 1 ) = i 

Thus the wave function at <£> = 2n does not change under the set of transformations 

( f) ( f) ( f) ( f) 

Pm ~* Pm + n an< ^ Pm-i ~^ Pm-i + 71 • Therefore the initial state and the final 

state are equivalent through the above set of the transformations. 
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APPENDIX C 

In this appendix, we show that the function (4.11) is odd function of $. We 
start from the twisted boundary condition 

M 

N Pj ($) + e(p j (*),p k (*))=2*I j + *, (C.l) 

fc=i(#) 

where N = 2M, Ij = 2j ~f -1 (j = 1, 2, • • • , M) and 

Asin[(j»j -pfe)/2] 



6(pj,Pk) = 2arctan 



cos[(pj+pk)/2] - Acos[(pj -pfe)/2] 



By noting 0(pj($), p fc ($)) = -0(p fc ($),Pj($)) and /, = -I M -j+l, (C.l)x(-l) is 
written as 

M 

N(- Pj (Q))+ 0(- Pj ^),-p k m = 2nI M+1 ^-^. 

Here we rewrite n + 1 — j to j; 

M 

Ni-pM+i-jW) + Yl 0(-PM+i-j($), -Pk(*)) = Stt/j - 

fe=l(^M+l-j) 
Compare this and boundary condition at (— $) 

M 

A%-(-ct>) + £ fl( Pi (-$),p Jfc (-$)) = 27r/ i -$, 
fc=i(#) 

we obtain 

Pj(-®) = -PM+i-#)- 
Using this relation and e(x^ — x ) = — e((— x a ) — (—#&)) the wave function (1.3) 
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becomes 



= Y\_ e ( X b ~ x a) ^(- 1) IQI e' ^ XaPQa ($)+ ^ ^°<fe & (P Q b ( $ )'^ (&))e(x b -x a ) 



a<b 



Q 



z (_ 1)M (M-l)/2 "Q e{{ _ Xb) _ { _ Xa)) ^ ( _ 1) |Ql e iE(-^)PM +1 - Qo (-*) 



x exp 



J2 d (PM+l-Q b (-$),PM+l-Q a (-®M(-Xb) - (~Xa)) 
a<b 



Here we put Q' a = M + 1 — Q a , then this wave function is expressed as 

= II e ((-^) - (-^)) ^(-l)IO'le i S(-^)Pg i (-*) e ^E.< fc «(PQi(-*)^(-*)) e ((-^)-(- 

a<6 Q' 

=x(-^i,---,-^m| -$), 

where we used (— l)'^' = (— 1)1*5 I x (— l) M ( M ~ 1 )/ 2 . Therefore we get the relation 

X(xi, . . .,xm\&) = x(~xi, -xm\ ~ 

From this relation and {xj} = {— xj} it is shown that the function (4.11) is odd 
function of $; 



E 
E 

{^} 



X(xi, . . . , x M \$)*S§x(xi, . . ■ , x M |$) 
||x(a:i...,XM|*)|| 2 

x(-xi, . . . , -x M | - $)*(-gf=$j)x(-zi, • • • > ~^m| - $) 
||x(-xi,...,-x M | -$)|| 2 



-E 



x(xi,...,x M | -$)*(g^y)x(a:i,---,a:M| - $) 
||x(xi,...,x M | -$)|| 2 
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